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Abstract 

We adapt the statistical mechanics of the shallow-water equations to the case where 
the flow is forced at small scales. We assume that the statistics of forcing is encoded 
in a prior potential vorticity distribution which replaces the specification of the 
Casimir constraints in the case of freely evolving flows. This determines a generalized 
entropy functional which is maximized by the coarse-grained PV field at statistical 
equilibrium. Relaxation equations towards equilibrium are derived which conserve 
the robust constraints (energy, mass and circulation) and increase the generalized 
entropy. 
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1 Introduction 



Two-dimensional turbulent flows with high Reynolds numbers have the strik- 
ing property of organizing spontaneously into large-scale coherent structures 
such as jets and vortices. Jovian atmosphere shows a wide diversity of struc- 
tures: Jupiter's great red spot, white ovals, brown barges,... A good hydrody- 
namical model to describe jovian atmosphere is provided by the Shallow- Water 
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(SW) equ ations. A statistical theory of the SW system has been developed re- 
cently by IChavanis fc Sommeria I (120021 ) . This extends the statisti c al me chan- 
ics of the incompress i ble 2D Euler equation proposed by iMiller I (Il990h and 



Robert fc Sommeria I (119911 ) to the case of compressible flows. It is however 



assumed that the flow is weakly compressible (small Mach number) so that 
the effect of waves is not dominant. Therefore, large-scale coherent vortices 
can form and persist for a long time. After a phase of chaotic mixing (violent 
relaxation), the system is expected to reach an equilibrium state which cor- 
responds to the most mixed state consistent with the constraints imposed by 
the dynamics. Mathematically, it is obtained by maximizing a mixing entropy 
while conserving energy, mass, circulation and all the higher moments of the 
PV vorticity (Casimirs). One difficulty with this approach is that the moments 
of PV depend on the resolution at which the PV field is considered so that 
the prediction can change accordingly. Furthermore, geophysical flows are usu- 
ally forced and dissipated at small-scales so that the cons ervation of all the 
Casimirs is abusive. In an attempt to solve these problems, lEUis et al I (120021 ) 



have proposed to fix a prior vorticity distribution instead of the Casimirs. It 
is assumed that this global distribution of vorticity is generated by the small- 
scale forcing so it must be given as an input in the statistical theory. In this 
approach, only the conservation of the robust constraints (E, T) is taken into 
account and the effect of the small-scale forcing is encapsulated in a prior 
vorticity distributi on or in a generaliz ed en tropy. This approa ch has been fur- 
ther developed in IChavanis I (j2005al ) and IChavanis I (j2005bl ) in the context 
of the quasi-geostrophic (QG) equations. In this paper, we show how it can 
be extended to the case of the Shallow- Water (SW) equations. Since this is 
only a slight, but interesting, var iant with respect to the un-fo rced case, we 
shall mostly refer to the study of IChavanis fc Sommeria I (120021 ) for technical 
details, without repeating all the steps of the derivation. 



2 The shallow-water equations 



The dynamical evolution of a thin fluid layer with thickness h(x, y, t) and 
velocity field u = (u,v)(x,y,t) submitted to a gravity acceleration g on a 
rotating planet is governed by the shallow water equations 

dh 

— + V-(/m) = 0, (1) 
du 

- + (u^V)u=-jV/i-2nxu. (2) 



The first equation can be viewed as an equation of continuity and the second 
equation as the Euler equation (in a rotating frame ft) for a barotropic fluid 
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with pressure p = \gh 2 . The Euler equation can be rewritten 
du 

-+(w + 2Q)xu=-V5, (3) 
at 

where w = ue z = V x u is the vorticity and where we have introduced the 
Bernouilli function 

B = ^ + gh. (4) 



The potential vorticity (PV) 

q = ""ft-' (5) 



is conserved for each fluid parcel, i.e. 

Each mass element /irfr is also conserved in the course of the evolution. This 
implies that the shallow-water equations conserve the PV moments 

r n = J q n hdv. (7) 



The moments n = 0,1,2 are, respectively, the total mass M, the circulation 
r and the PV enstrophy T 2 . The energy 

E = J h^dr gh 2 dr, (8) 



involving a kinetic and a potential part is also conserved. It is convenient to 
use a Helmholtz decomposition of the momentum hu into a purely rotational 
and purely divergent part 

hu = -e z x + V0. (9) 



For any stationary solution, the mass conservation (1) reduces to V • (hu) = 
so that 

hu = -e z x W, (10) 
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where if) is the stream-function. Then Eq. (6) reduces to u • Vq = which 
implies that q = f(ip). Finally, Eq. (3) reduces to 



(u + 2J])xu = -VB. (11) 



Taking the scalar product with u, we obtain u • V5 = so that B = B(if>). 
Then, substituting Eq. (10) in Eq. (11) we obtain 

dB . . 

' = -*■ (12) 



The SW equations admit an infinite class of stationary solutions, specified by 
the relations B = B(if>) and q = f{if>) = —B'{ifi). They are determined by the 
two coupled partial differential differential equations for ip and h 

h 2 h h 3 dip ' 



2h 2 



+ gh = B(i/j). (14) 



This formulat ion of the SW equations in t erms of (h, if}) variables has been 



introduced in IChavanis fc Sommeria I (120021 ). 



3 The equilibrium statistical mechanics 

3. 1 Freely evolving flows 



The SW equations are known to develop a complicated mixing process which 
ultimately leads to the emergence of a large-scale coherent structure, typically 
a jet or a vortex. One question of fundamental interest is to understand and 
predict the structure and the stability of these equilibrium states. This can 
be achieved by using statistical mechanics arguments. The idea is to replace 
the deterministic description of the flow q(r, t) by a probabilistic description 
where p(r, a, t) gives the density probability of finding the PV level q = a 
in r at time t (it satisfies the local normalization condition / p(r,o~)do~ = 1). 
The observed (coarse-grained) PV field is then expressed as q(r, t) = J pada. 
To apply the statistical theory, one must first specify the constraints attached 
to the SW equations. The mass M = J hdr, circulation r = / qhdr and 
energy E = J h^-dr + k I gh 2 dr will be called robust constraints because they 
can be expressed in terms of the coarse-grained field. These integrals can be 
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calculated at any time from the coarse-grained field and they are conserved 
by the (macroscopic) dynamics. By contrast, the Casimir invariants If = 
I f(l)hdr, or equivalently the fine-grained moments of the vorticity = 
Jq^hdr = J pa n dahdr, will be called fragile constraints because they must 
be expressed in terms of the fine-grained PV. Indeed, the moments of the 
coarse-grained PV r^>i = / cphdr are not conserved since 7^ q n (part of 
the coarse-grained moments goes into fine-grained fluctuations). Therefore, 
the moments r{>i must be calculated from the fine-grained field q(r, t) or 
from the initial conditions, i.e. before the PV vorticity has mixed. Since we 
often do not know the initial conditions nor the fine-grained field, the Casimir 



invariants often appear as "hidden constraints" (IChavanis 1 . 120061 ). 



The statistical theory of Miller and Robert-Sommeria-Chavanis for the 2D 
Euler equations and the SW equations is based on several assumptions: (i) 
it is assumed that the flow is freely evolving without small-scale forcing and 
dissipation, (ii) it is assumed that we know the initial conditions (or equiva- 
lently the value of all the Casimirs) in detail, (iii) it is assumed that mixing is 
efficient and that the evolution is ergodic so that the system will reach at equi- 
librium the most probable (most mixed) state. Within these assumptions^.], 
the statistical equilibrium state of the SW system is obtained by maximizing 
the mixing entropy 

S[p] = — / pin phdrda, (15) 



while conserving the energy, the circulation, the mass and all the Casimirs. 
We write the variational principle in the form 

5S - f35E - a5T - 7 5M - £ aJT f n 9 - - J ((r)5 (J pdaj hdr = 0. (16) 



In the present point of view, all the constraints are treated on the same footing. 



1 Some attempts have been propo sed to go beyond the assumptio ns (ii) and (iii) of 
the statistical theory. For example, Chavanis Sommeria ( 19961 ) consider a strong 
mixing limit in which only the first moments of th e vorticity are relevant i nstead 
of the whole set of Casimirs. On the other hand, IChavanis Sommeria ] (|l998h 
introduce the concept of maximum entropy bubbles (or restricted equilibrium states) 
in order to account for situations where the evolution of the flow is not ergodic 
in the whole available domain but only in subdomains. A 2D turbulent flow is 
therefore viewed as an ensemble of isolated vortices which can be seen as "maximum 
entropy bubbles" separated by un-mixed flow. In 2D decaying turbulence, these 
isolated vortices result from prev ious mergings and they are expected to correspond 



to statistical equilibrium states ( Laval et al . 200ll ). In fact, because of incomplete 



relaxation, they may just well be particular stable stationary solutions of the 2D 
Euler equation that are incompletely mixed. 
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In particular, the moments T^ 9 ' are treated microcanonically and we must 
ultimately relate the Lagrange multipliers a n to the constraints r£ s \ 



3.2 Prior distribution and relative mixing entropy 



In the approach of IChavanis fc Sommeria I (120021 ). it is assumed that the sys- 
tem is strictly described by the SW equations so that the conservation of all the 
Casimirs has to be taken into account. However, in geophysical situations, the 
flows are forced and dissipated at small scales (due to convection i n the jovian 



atmos phere) so that the conservation of the Casimirs is destroyed. lEUis et al 



( 120021 ) have proposed to treat these situations by fixing the conjugate vari- 
ables a n> i instead of the fragile moments If we view the PV levels as 
species of particles, this is similar to fixing the chemical potentials instead of 
the total number of particles in each species. Therefore, the idea is to treat the 
fragile constraints canonically, whereas the robust constraints (E, T, M) are 
still treated micro cano nically. Thi s point of view has been further developed in 
Chavanis I (j2005al ) and IChavanis I (j2005bl ) in the QG model and this approach 
is extended here to the SW system. The relevant thermodynamical potential 
is obtained from the mixing entropy (15) by using a Legendre transform with 
respect to the fragile constraints: 



5 V = S 



(17) 



n>l 



Explicating the fine-grained moments, we obtain 



S x = - p In p + ctncr n 



hdrda. 



Introducing the function 
X((r) = expi - J2 a n<? T 

n>l 



(19) 



we get 

s x [p] 



p In 



hdrda, 



(20) 



which has the form of a relative entropy. The function x{ a ) is interpreted as a 
prior vorticity distribution. It is a global distribution of PV fixed by the small- 
scale forcing. We shall assume that this function is imposed by the small-scale 
forcing so it must be regarded as given. 
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Assuming ergodicity, the statistical equilibrium state is now obtained by maxi- 
mizing the relative entropy S x at fixed energy E, circulation T and mass M (ro- 
bust constraints). The conservation of the Casimirs has been replaced by the 
sp ecification of the prior x ( <?)- Fr om that point, we can repeat the calculations 
of Chavanis fc Sommeria I ( 2002 ) with almost no modification. The only dif- 



ference is that we regard the a n as given. Writing 5S X — (35 E — a5T — ^bM = 0, 
and accounting for the normalization condition / pdo = 1, we find (i) that the 
statistical equilibrium state is a stationary solution of the SW equation (ii) 
that the detailed distribution of PV levels is given by the Gibbs state 



P(r,<r) 



1 



-{(3ip+a)cr 



(21) 



where the partition function 



(22) 



is determined by the local normalization condition. The distribution of the 
fluctuations of PV vorticity is the product of a universal Boltzmann factor by 
a no n-universal function y{cr) fixed by the forcing. This is the same formula 



as m 



Chavanis fc Sommeria I (120021 ) except tha t in the present formalism y(cr 



must be regarded as given a priori while in IChavanis fc Sommeria I (120021 ) 
it was a function of the Lagrange multipliers a n that had to be related a 
posteriori to the fine-grained moments imposed by the initial conditions. The 
equilibrium coarse-grained PV is 



/ xi^o-e-^+^da 1 din Z 
fx(o-)e-(M+ a ^da ~ ~~(3 di\) 



*W + «) = /(# 



(23) 



The coarse - grain ed vorticity (23) can be viewed as a sort of super statistics 



(IChavanis 1 . 120061 ) as it is expressed as a superposition of Boltzmann factors 
(on the fine-grained scale) weighted by a non- universal function We note 
that the q — ip relationship predicted by the statistical theory can take a wide 
diversity of forms (usually non-Boltzmannian) depending on the prior x{ a )- 
The function F is entirely specified by the prior PV distribution according to 



= -Qn *)'($), with = J X (v)e- a *da. 



(24) 



Differentiating Eq. (23) with respect to ip, we find that 



q'(ip) = -f3q 2 , q2 = q 2 ~q 2 > 0. 



(25) 
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This relation relates the slope of the q = f{ip) relationship to the local 
centered variance q^ (ip) of the PV distribution. This can also be written 
F'(<&) = — q2^) < so that F is a decreasing function. Therefore, the sta- 
tistical theory predicts that the coarse-grained field is a stationary solution 
of the SW equation and that the q — ip relationship is a monotonic function 
which is increasing at negative temperatures (3 < and decreasing at posi- 
tive temperatures (3 > 0. We note that, according to Eqs. (12) and (23), the 
Bernouilli function is given by 

B = ^\nZ, (26) 

and it plays the role of a free energy in the statistical theory (if we interpret 
Z as a partition function). 

We also note tha t the most probable PV field (a) (r) of the Gibbs distribution 



(21) is given by (ILeprovost et al 1 . 120061 ): 



(a) = [(ln X )rW + «), (27) 

provided that (lnx)"((cx)) < 0. This is also a stationary solution of the SW 
system which usually differs from the average value q(r) of the Gibbs distri- 
bution (21). They coincide only when 

-(lnx)'($) = [(ln X )r 1 ($), (28) 
which is the case when x(°") is gaussian. 

3.3 Generalized entropy 

We first show that a PV field which extremizes a functional of the form 

H[q] = - I C(q)hdr, (29) 



where C is a convex function, at fixed energy E, mass M and circulation T, 
is a steady state of the SW equations. We write the variational principle as 

5H - I35E - a5T - j8M = 0. (30) 



Using the results of IChavanis fc Sommeria I (120021 ). we have 



5E = J B5hdv + J ip5(qh)dr - J 05(V • u)dr, (31) 
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ST = J 5(qh)dr, (32) 
SH=- J C{q)5hdr - J C'(q)8(qh)dr + J C'{q)q5hdr, (33) 



where we have taken h, qh and V • u as independent variables. The variations 
on V • u yield = 0. The variations on qh give 

C\q) = -(3ip - a, (34) 



so that q = f(tp). The variations on h give 

qC'(q)-C(q)=(3B + 7} (35) 



so that B = B{ip). Taking the derivative of Eq. (35) with respect to ip, we 
find that 

" C "^Ti = C- < 36 > 

dip dip 



According to Eq. (34), we also have 

dq_ 
# ~ ~C"{qY 



(37) 



so that q(ip) is a monotonic function increasing at negative temperatures and 
decreasing at positive temperatures. Furthermore, we find that Eq. (36) is 
equivalent to q — —dBjdip. Therefore, the optimization problem (30) de- 
termines stationary solutions of the SW system. Since H, E, T and M are 
conserved by the SW equations, we can argue, as for the 2D Euler equation 



(jEllis et al 1 . 12002| ). that a maximum of H at fixed E, T and M (if it exists) 
will be nonlinearly dynamically stable with respect to the SW system. In this 
dynamical context, H is referred to as a Casimir functional and E — H as an 
energy-Casimir functional. 

Note that the optimization problem (30) can also be justified by a selective 
decay principle (for —H) if we interpret the PV vorticity as the coarse-grained 
PV. Indeed, —H[q] calculated with the coarse-grained PV decreases (fragile 
constraint) while E\g], T[q] and M\q] are approximately conserved (robust 
constraints). This selective decay principle can explain physically how —H[q] 
can possibly reach a minimum value while —H[q\ is exactly conserved by 
the SW equations. In t his coarse-grained con text, H\q] is referred to as a 



generalized if-function (jTremaine et al I . [1986J ). 
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Finally, we note that the equilibrium state (23) predicted by the statistical 
theory extremizes a certain if -function at fixed E, F and M. This functional 



S[q] = -jc(q)hdr, 



(3* 



corresponding to t he statistical equi librium state, will be called a generalized 



entropy in g-space (IChavanis 1 . 120061 ). It is completely determined by the prior 



PV distribution. It should not be confused with the mixing entropy (15) which 
is a functional of p. Coming back to Eq. (34) and recalling that C is convex (so 
that this relation can be inversed), we find that Eqs. (34) and (23) coincide 
provided that C'(q) = — F _1 (g). Therefore, the prior PV distribution x(cr) 
determines F(x) which itself determines C(q) according to 



C(q) = - j F-\x)dx. 



(39) 



In other words, to obtain C(q), we need to inverse Eq. (23) and integrat e 
the resulting expression. Some examples are collected in IChavanis I (120031 ). 
Combining the previous relations, we find that the g eneralized entropy is de- 
termined by the prior according to (IChavanis 1 . 120061 ) 



C{q) 



[(Inx) 



'i-i/ 



-x)dx. 



(40) 



Finally, comparing Eqs. (37) and (25) we get the relation 
1 



Q'2 



C"{q) 



(41) 



which is exact at statistical equilibrium. 



4 Relaxation equations 



4-1 Maximum Entropy Production Principle 



In the case of freely evolving flows, IChavanis fc Sommeria I (120021 ) have pro- 
posed a thermodynamical parametrization of the SW equations (on the coarse- 
grained scale) in the form of relaxation equations that conserve all the con- 
straints of the SW dynamics (including the Casimirs) and increase the mixing 
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entropy. In the case of flows that are forced at small-scales, the Casimirs are re- 
placed by the specification of a prior vorticity distribution or, as we have seen, 
by a generalized entropy. In this context, we can propose a thermodynami- 
cal parametrization of the SW equations in the form of relaxation equations 
that conserve only the robust constraints (mass, energy and circulation) and 
increase the generalized entropy S[q] fixed by the prior. 

We first decompose the vorticity u and velocity u into a mean and fluctuating 
part, namely uj = u + u, u = u + u, keeping h smooth. Taking the local 
average of the shallow water equations (1)(3), we get 

dh 

- + V • (hu) = 0, (42) 

- + (u+2fi)xu=-VB-e z xJ U) (43) 



where the current J w = cDu represents the correlations of the fine-grained 
fluctuations. We deduce an equation for the evolution of the potential vorticity 
(5), taking the curl of Eq. (43) and using Eq. (42): 

^(/ig)+V-(/igu) = -V.J w . (44) 



This equation can be viewed as a local conservation law for the circulation 
r = / qhdr. It shows also that represents the current of coarse-grained 
vorticity due to mixing. We shall determine the unknown current by the 
thermodynamic principle of Maximum Entropy Production (MEP) , using the 
generalized entropy (38). The Maximum Entropy Production (MEP) princi- 
ple consists in choosing the current which maximizes the rate of entropy 
production S respecting the constraints E — 0, and < C(r,t). The last 
constraint expresses a bound (unknown) on the value of the diffusion current. 
Convexity arguments justify that this bound is always reached so that the 
inequality can be replaced by an equality. We write the variational problem 
as 

5S - p(t)6E- J D- 1 5(^ S )dr = 0, (45) 



where f3(t) is a Lagrange multiplier accounting for the conservation of energy 
and D~ l is a Lagrange multiplier associated with the constraint = C(r,t). 
The conservations of mass and circulation are automatically satisfied by the 
form of the relaxation equations (42) and (43). Noting that 

E = fj u - u ± hdr, (46) 
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S = - J C"(q) J. • Vgrfr, 



(47) 



and performing the variations on in Eq. (45), we obtain an optimal current 



-D 



(4£ 



Thus, in the presence of a prior PV distribution, we obtain a parametrization 
of the SW equations of the form 



^ + V ■ (hu) = 0, 
^- + qhe z xu=-v(y + 3/i) +D 



e z x Vg 



C"(9) 



hu 



_ (V x u) ■ e z + 2n 

Q = 7 , p(t) 



J Dhu ± ■ Vqdr 



n • ui (on each boundary), 
C"(g) 7 



n ■ Vg = 

n • u = (on each boundary) 



(49) 
(50) 

(51) 

(52) 
(53) 



where n is a unit vector normal to the boundary and we have omitted the 
overbar on u. The relaxation equation for the coarse-grained vorticity is given 
by 



^(hq)+V-(hqu)=V-lD 



Vq 



Q 



— huj_ 



(54) 



These equations can al s o be directly obtained from the parametrization of 
Chavanis &: Sommeria I (120021 ) by replacing q2 in their parametrization by 
1/C"(q). Therefore, the identity (41) can be viewed as a closure relationship 
in the present context. This relation is valid at equilibrium but the present 
approach suggests that it is also valid out-of-equilibrium when there is a prior 
distribution of PV. In fact, we can obtain this relation by assuming that, 
out-of-equilibrium, the PV distribution p(r, a, t) maximizes the r elative en 



tropy ( 20) at fixed PV q(r, t) and normalization (Appendix C of I Chavanis 
(l2005i 



The entropy production (47) can be written 



S = - C"(q)3 u 



dr + (3(t) / J w ■ u±hdr. 



(55) 



12 



Using the conservation of energy E = with Eq. (46), the second integral is 
seen to vanish. Substituting Eq. (48) in the first integral, we finally obtain 

S = J C"(q)^dr>0, (56) 



which is positive provided that D > 0. A stationary solution of the relaxation 
equations (49)- (53) satisfies S = yielding = 0, i.e. 

w+ c4) w=o ' (57) 



After integration, we obtain 

C'(q) = -I3tfj - a. (5S 



Therefore, the generalized entropy increases until the statistical equilibrium 
state (23)-(34), fixed by the prior xO 7 )) is reached with (3 = lim^ +00 (3{t). 
Alternatively, these equations can be used as a numerical algorithm to com- 
pute arbitrary stationary solutions of the SW system specified by the convex 
function C (see Sec. 3.3). 



4-2 The incompressible limit 



The case of ordinary 2D incompressible turbulence is recovered in the limit 
h — > 1, q — > u) and u = — e z x V^. The relaxation equation for the coarse- 
grained vorticity is given by 



duJ „ 



C"{TD) 



(59) 



with 



(3{t) 



J DVuJ ■ Vtpdr 

J U C"(TJ) al 



(60) 



This returns the equations introduced by IChavanis I (120031 ) in 2D turbulence 
in the case where the system is described by a prior vorticity distribution. 
They can be viewed as nonlinear mean-field Fokker-Planck equations. They 
are the forced-case counter part of the relaxation equations introduced by 



Robert fc Sommeria I (119921 ) for freely evolving flows that conserve all the 



Casimir constraints. 
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The relaxation equation (50) for the velocity field can be written 



c)w 1 

-t^ + (u- V)u= — Vp + D 

at p 



A W 



(61) 



where we have used the well-known identity of vector analysis Au = V(V • 
u) — V x (V x u) which reduces for a two-dimensional incompressible flow to 
Au = e z x Vu7. Eq. (61) is valid even if D is space dependant unlike with a 
usual viscosity term. In previous publications this equation was given only in 
its vorticity form (59) and the equivalence with Eq. (61) is not obvious at first 
sights when D is space dependent. At equilibrium, we have from Eq. (61) the 
identity 



Au 



u. 



C"(uJ) ' 



(62) 



which can be deduced directly from the stationary state (34). Indeed, for a 
stationary solution To = f(ip), the previous identity Au = e z x Vw becomes 
Au = —f'(ijj)u which is equivalent to Eq. (62) for a steady state thanks to Eq. 
(37). Finally, the previous equations can be extended to the quasi-geostrophic 
(QG) limit if we replace uj by the PV q related to the stream-function by 

q = -A^ + (63) 



where L R is the Rossby radius. 



4-3 Explicit examples 



In the present formalism, one has to specify (i) a prior PV distribution x(cr) 
which encodes small-scale forcing and non-ideal effects (ii) the robust con- 
straints E, T and M which can be determined at any time from the coarse- 
grained flow. From the prior PV distribution, we can determine a generalized 
entropy C(q) by using the procedure exposed in Sec. 3.3 (see in particular 
Eq. (40)). Then, we can use this entropy in the parametrization (49)- (53) to 
obtain the dynamical evolution of the coarse-grained flow towards statistical 
equilibrium. The difficulty is now to find the good prior. This depends from 
case to it is related to the properties of forcing, but the idea is that 

several forms of prior (or corresponding entropies ) give simi l ar res ults so that 



they can be regrouped in classes of equivalence (IChavanis 1 . 120031 ) . Thus, for 



a given situation, one has to find the relevant class of equivalence. In general, 
one has to proceed by tryings and errors. We specify a prior, compute the 
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corresponding flow and see whether it agrees with the information that we 
have on the system. If we find a "good prior" for a system with some given 
initial conditions, then we can expect that it will remain valid for this system 
when we change the initial conditions (or, equivalently, the value of the robust 
constraints E, T and M). Note that specifying the prior is not the end of the 
story but only the starting point. Indeed, there can be different types of solu- 
tions for a given prior PV distribution depending on the values of the control 
parameters. Thus, we have to study the bifurcation diagram as a function of 
these parameters E, F and M for a given prior x( cr ) or generalized entropy 
C(q). This is a rich and non trivial problem. 



Let u s give some examp les to illustrate our approach. In the case of Jovian 



flows, Ellis et al I (120021 ) have proposed to adopt a prior x( a ) corresponding 



to a de-ce ntered Gamma d istribution. This leads to a generalized entropy of 



the form (IChavanis 1 . 120061 ) 



C(q) = ~ 



q ln(l + eg) 



(64) 



where 2e is the skewness of the PV distribution. We have C'(q) = q/(l + eg) 
and C"(q) = 1/(1 + eg) 2 so that the statistical equilibrium state is specified 
by 

(65) 



In the limit e — > 0, the generalized entropy (64) becomes minus the enstrophy 
S[q] = — \ J q 2 dr and the q = f(ip) relationship is linear. The parametrization 
that we propose in that case for the vorticity equation is 



d 

— (hq)+V-(hqu) 



v- Id 



Vq + (3(t)(l + eqYhu ± 



(66) 



On the other hand, in order to describe jovian flows, ISommeria et al I (Il99lh 



have considered of the statistical theory where the PV distribution is 

restricted to two-levels. In their approach, the flow is assumed to be freely 
evolving and the dynamics corr esponds to t he invi scid mixing of patches with 
P V cr an d o"! . As discussed in IChavanis I (j2005a ) , we can also consider the 
case of a flow that is permanently forced at small-scales (due to convection) 
so that a prior PV distribution is established with two intense peaks at cr and 
a\. In the two-levels case, theses two interpretations lead to the same results 
but the second one can probably be extended more easily to the more realistic 
situation where the peaks have a finite width Act. If we adopt a prior of the 
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form x( a ) = fi( a ~~ + X^i 17 ~ a o), we find that the generalized entropy is 



C(q) 



a i - cr Q 



[(q - (T ) \n(q - <x ) + - q) ln(<Ti - g)] 



(67) 



We have C"(g) 

statistical equilibrium state is 



^— ln(^) and C"( 

<xi— o"o CT i — Q 



l/[(g — cto)(o"i — g)] so that the 



q = a + 



(6S 



similar to the Fermi-Dirac distribution. The parametrization that we propose 
in that case for the vorticity equation is 



Of 



(hq) + V ■ (hqu) = V -\D 



Vq + (3{t){q - a )(a 1 -q)hu ± 



(69) 



and it coincides with the parametrization of IChavanis fc Sommeria I (120021 ) 
in the two-levels case of the statistical theory Therefore, the approach based 
on priors is not in radical opposition with the usual statistical theory but it 
allows for convenient extensions in more general cases. 

Explicit determinations of the statistical equilibrium state specified by Eqs. 
(65) and (68) have been obtained in the QG limit of the statistical theory. 
In the two approaches, a vortex solution has been found at the latitude of 
Jupiter's great red spot where the underlying topography is extremum (or 
the shear is equ a l to z ero). In the case of the Fermi-Dirac distribution (68), 
Sommeria et al I (Il99ll ) consider the limit of small Rossby radius and find 
that the vortex has an annular jet structure which is consistent with the 
m orphology of Jupiter's g r eat re d spot. This study has been fur ther developed 



in 



Bouchet fc Sommeria I (120021 ) and it has been extended by iBouchet et al 



(120061 ) to the case of the SW system. In these studies, the structure of Jupiter's 
great red spot can be seen as the co-existence of tw o thermodynamica l phases 



in contact separated by a sort of "domain wall" (IChavanis I . l2005al ) in the 



langage of phase ordering kinetics (like in the Cah n-Hilliard theory) . This 
annular jet structure is not obtained in the approach of lEUis et al I (120021 ). This 
implies that the prior PV distribution relevant to the case of JGRS is probably 
closer to two intense peaks rather than to a decentered Gamma distribution. 
These two distributions belong to different classes of equivalence since they 
lead to structurally different solutions. Probably, the prior distribution could 
be improved to give a better description of jovian vortices, but a distribution 
with two intense peaks already gives a fair description. 



16 



5 Conclusion 



In this paper, we have extended the "ordinary" statistical theory of the SW 



system (IChavanis fc Sommeria 1 . 120021 ) to account for the existence of a prior 



vorticity distribution. This approach can be justified when the system is forced 
at small-scales so that a permanent PV distribution is imposed canonically. 
In that case, the forcing acts like a sort of reservoir. The attractive nature of 
this theory is its practical interest: in the standard theory, one works with the 
Boltzmann entropy in p-space and deals with a very large (possibly infinite) 
number of Casimir constraints which are often not known or not rigorously 
conserved. In the other approach, one conserves only the robust constraints 
(E, T, M) and work with a generalized entropic functional in g-space fixed by 
the prior (that has to be found by tryings and errors). Therefore, in the first 
approach, we have to solve N coupled relaxation equations (one for each level) 
while in the second approach, we just have to solve one relaxation equation. 
Whether this approach is really physically relevant remains to be established. 
In any case, the relaxation equations (49)-(53) can be used as numerical algo- 
rithms to construct a large class of stationary solutions of the SW equations, 
which is certainly an interest of our formalism. 



References 

Bouchet, F., Sommeria, J.: 2002, JFM 464, 165 

Bouchet, F., Chavanis P.H. & Sommeria, J.: 2006, in preparation 

Chavanis, P.H.: 2003, PRE 68, 036108 

Chavanis, PH.: 2005a, Physica D 200, 257 

Chavanis, P.H.: 2005b, in Proceedings of the Workshop on "Interdisciplinary 
Aspects of Turbulence" at Ringberg Castle, Tegernsee, Germany (Max- 
Planck Institut fiir Astrophysik) (physics/0601087 

Chavanis, PH.: 2006, Physica A 359, 177 

Chavanis, PH., Sommeria, J.: 1996, JFM 314, 267 

Chavanis, PH., Sommeria, J.: 1998, JFM 356, 259 

Chavanis, PH., Sommeria, J.: 2002, PRE 65, 026302 

Ellis, R., Haven, K., Turkington, B.: 2002, Nonlinearity 15, 239 

Laval, J.P, Chavanis, PH., Dubrulle, B., Sire, C: 2001, PRE 63, 065301(R) 

Leprovost, N., Dubrulle, B., Chavanis, PH.: 2006, PRE 73, 046308 

Miller J.: 1990, PRL 65, 2137 

Robert, R., Sommeria, J.: 1991, JFM 229, 291 

Robert, R., Sommeria, J.: 1992, PRL 69, 2776 

Sommeria, J., Nore, C, Dumont, T., Robert, R.: 1991, CRAS Ser. II 312, 999 
Tremaine, S., Henon, M., Lynden-Bell, D.: 1986, MNRAS 219, 285 



17 



